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O Abstract 

Renormalization group limit cycles may be a commonplace for quantum 
Hamiltonians requiring renormalization, in contrast to experience to date with 
classical models of critical points, where fixed points are far more common. 
We discuss the simplest model Hamiltonian identified to date that exhibits a 
renormalization group limit cycle. The model is a discrete Hamiltonian with 
two coupling constants and a non-perturbative renormalization group that in- 
volves changes in only one of these couplings and is soluble analytically. The 
^ Hamiltonian is the discrete analog to a continuum Hamiltonian previously 

^ , proposed by us. 



1 Introduction 

In 1971, one of us suggested that renormalization group equations could have limit cycle 
solutions as well as fixed points when the renormalization group equations involve two 
or more coupling constants [l|]. The 1971 paper did not mention the possibility that 
renormalization group equations might have limit cycle solutions even for differential 
equations for only one coupling constant. But in 1993, the two of us defined a simple 
Hamiltonian that requires the renormalization of just a single bare coupling constant <?a, 
where A is the cutoff @. The bare coupling constant, we demonstrated, exhibits limit 
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cycle behavior in the limit A — > oo. To be precise the coupling constant g\ was found 
to decrease steadily as A increased, until g\ reached — oo, after which it jumped to +00 
and started its next cycle of steady decrease. However, we did not recognize or comment 
that this behavior constituted a limit cycle. Moreover, the Hamiltonian we defined is 
singular, in such a way that it presumably does not have well-defined s-wave scattering 
phase shifts. The model has well-defined bound states, and our renormalization group 
analysis was based on keeping the finite bound state energies fixed as A — > 00. 

The model Hamiltonian we published in 1993 has no known physical application. 
But in 1999, Bedaque, Hammer, and Van Kolck showed that a three-body Hamiltonian 
with two- and three-body delta function potentials and a cutoff in momentum space is 
renormalizable and that the three-body coupling approaches a limit cycle as the cutoff 
lambda approaches 00 [||], much as it does in our 1993 model. Bedaque et al built on 
earlier work of Thomas |3| and Efimov || , that was recently reviewed by Nielsen, Fedorov, 
Jensen, and Garrido ||. However, the mathematics of the three body system is complex, 
and it is surely useful to have simpler models to study that exhibit limit cycles too. 

In this paper, we discuss a discretized version of the Hamiltonian we introduced in 
1993. With a cutoff, the discretized model takes the form of a finite size matrix with 
discrete eigenvalues. The discrete matrix has a diagonal sub-matrix, plus two off-diagonal 
pieces with two coupling constants. The Hamiltonian requires renormalization in the limit 
of infinite cutoff. Just as in the continuum case, the renormalization can be constructed 
analytically and leads to a limit cycle in one of the two couplings, while the second 
coupling constant stays fixed. But the analytically obtained limit cycle is defined only 
for a discrete sequence of cutoffs, rather than for a continuously varying cutoff. Using 
numerical procedures, it is possible to define and compute a renormalization process with 
a continuously varying cutoff, but our studies with the continuous cutoff variation will 
not be discussed here. 

The finite matrix Hamiltonians can be diagonalized numerically, when the cutoff is 
small enough. The renormalizability of the Hamiltonians we will discuss here can be 
demonstrated to high numerical accuracy with cutoffs small enough to allow numerical 
diagonalization. We will provide the demonstration with a comparison of the eigenvalues 
of two matrices with two different cutoffs, one of size 37x37, the other of size 42x42. 

The Hamiltonian to be used in this paper has the form 



where m and n are integers. For m = n, 5 mn = 1 and s mn = 0. For m 7^ n, 5 mn = and 
s mn = (m — n)/\m — n\. The numbers E n = b n with b > 1, are eigenvalues of the operator 
Hq that has matrix elements (m|f/o|^) — H mn (0,0). The eigenvalues are called kinetic 
energies of the corresponding eigenstates, |n). The remaining part of the Hamiltonian, 
Hj = H — H , is called an interaction, and Hj(0, 0) = 0. The largest energy allowed 
in the dynamics is A N = b N , which defines the ultraviolet cutoff so that the subscripts 
m, n < N. 

The continuous version of this model is recovered in the limit b —*■ 1. The discrete 
model itself has been discussed in the case = using similarity renormalization group 
idea M and Wegner's equation ||. It was shown that the model exhibits asymptotic 
freedom for = 0. Hamiltonians of Eq. ([!]) with = can be derived in a number of 



H m n(gN, h^) — (E m E n ) 
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ways, ranging from a discretization of a nonrelativistic Schrodinger equation for a particle 
on a plane with a two-dimensional ^-potential, to a discrete version of the transverse 
dynamics of partons in quantum field theory. 

At first, the model with ^ does not appear much different from the one with 
hpf = 0. All Hamiltonians defined by Eq. (|1|) are hermitian and have a general ul- 
traviolet logarithmically divergent structure. The divergence originates in the fact that 
the far off-energy- shell matrix elements of H mn , i.e. those with \m — n\ ^> 1, behave 
like (EmEn) 1 / 2 . Therefore, perturbation theory produces large-energy contributions of 
the form Y.i(EmEi) 1 ^ 2 {1/Ei) (EiEn) 1 / 2 , and every energy scale contributes a significant 
amount to the sum, which grows linearly with the number of initially included scales, 
iV = In Aw/ In 6. The dependence of the unrenormalized coupling constants and h N 
on N, is expected to compensate this effect. As we will demonstrate below, when is 
not zero, the model exhibits limit cycle behavior as N goes to infinity. We will prove this 
by deriving a renormalization group equation that determines Qn-x (used with a cutoff 
Ajv-i ), given g N and h N used with cutoff A N , such that the low energy eigenvalues stay 
fixed. 

This paper is organized as follows. Section describes the discrete RG flow of the 
model Hamiltonians H(gN,htf) with N. The structure of eigenstates of these Hamilto- 
nians is also briefly described. Section || shows a numerical example of the spectrum of 
renormalized Hamiltonians and their discrete rescaling properties. The spectrum is found 
to contain a sequence of negative (bound state) eigenvalues, one per cycle. Section | 
concludes the paper. 

2 Renormalization group trajectories 

The eigenvalue problem 

N 

H mn tp n = Eip m , (2) 

n=— oo 

can be solved for ip m , m < N, assuming that one knows E, by using the Gaussian elimi- 
nation procedure. In the first step one can express ^ in terms of all other components, 
ip n with n < N. In the next step, one expresses ipN-i m terms of components ip n with 
n < N — 1, and so on. In fact, this procedure can be carried out without knowing E if 
the eliminated energy scales are much larger than E. Thus, when one is interested in the 
eigenvalues of order b M with M — > — oo, one can eliminate states from some huge energy 
range between = b N and some fixed A = b N °, where N is much smaller than N and 
still A 3> E. The new, much smaller eigenvalue problem will have the cutoff A and 
the corresponding coupling constants. When the recursion for the whole sequence of the 
coupling constants (#jv, /iat), (pjv-i, ^at-i),--, (9N ,hN ), is found, one can trace it back 
starting from some assumed finite values (go, ho) at a finite A^ and see what happens with 
(gN, h N ) when the cutoff number — > oo. 

To carry out the above procedure, it is convenient to introduce a small number E = 
E/b N and write the eigenstate components in the form ip n = b n l 2 <\) n for all n < N. The 
eigenvalue Eq. (^) becomes: 
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(1 - E)<p N = g N a N + ih N a N ^i , 
(1-6^)0jv-i = 9nOn + ih N {-<I>n + crjv-2) , 

1 - b 2 E)(f) N ^ 2 = 9N<?N +ih N (-(j) N - 0jv_i + fXAr_ 3 ) 



1 - b n E)<p N _ n = g N a N + ih N (-<p N + — - <pN-n+i + ^N-n- 



1) > 



where <7j = S] = _ DO 0j. If and n are finite and A" — > 00, fe*^ for i = 0, 1, n on the left 
hand sides of the above equations do not count. The first step of the Gaussian procedure 
gives 

<Pn = (gN + ih N )a N -t/(l - g N ) . (3) 

Substituting this result into the remaining equations, one obtains a new set that does not 
explicitly involve the component with largest kinetic energy but has a different coupling 
constant instead. Namely, 

(1 - bE)(p N ^i = gN-i^N-i + ih N _i(T N _ 2 , 

(1 -b 2 E)(j) N ^2 = gN-X°N-l + itlN-l (-0AT-1 + CrjV-3) , 

(1 -b n E)(j) N _ n = g N -l(?N-l +ih N -l {~<f>N-l + ■■■ ~ ^N-n+l + CjV-n-l) , 



The new coupling constant gN-i = 9n + (g% + h N )/(l — g N ), while h^-i = h N stays 
unchanged. Making p such steps, one transforms (gN,h^) into (g^- p ,hisr- p ), where 
h-N-p = and 



9n- p = h N tan 



arctan ( |— ) +p arctan(Zijv) 
h N J 



(4) 



It is visible that g^- P = gN when = tan(n/p). Equation (§) demonstrates that the 
coupling constant g^ completes a cycle when A^r is reduced to A^/lf, and the cycle repeats 
itself in the further reductions. Fig. 1 shows these cycles for Ai 6 = b 16 , g = 0.060606 |§, 
and h w = tan(7r/5). If g±e were zero, the cycle would be symmetric with respect to the 
axis gpf = 0. On the other hand, components of the eigenvectors of H(g N , h^) satisfy the 
recursion 

1/2 l-ih N -E/b n+1 

^n = b> — — i> n+1 , 5 

1 + ih N — E jb n 

for all A^ and n < N independently of the value of g^- 

Note that for real Hamiltonian matrices, i.e. for hjsi — > 0, Eq. (^implies 

1 +g (N - N ) 

which means asymptotic freedom, or g& = g / [1 + g ln(A/A )/ln6] — > when A — > oo. 
This is a fixed point behavior. In contrast, for complex Hamiltonians, Eq. (P is equivalent 
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Figure 1: The discrete limit cycle dependence of the coupling constant g n on the cutoff 
number iV = lnA^/lnfe, g 16 = 0.060606, h\ 6 = tan(7r/5). The lines are added between 
the discrete points only to guide the eye. 

to the following renormalization group transformation of g\ from A to A, (h = = 
const.), 

_ , go ~ h tan [cln(A/A )] 
9A ~ h + g tan[cln(A/A Q )j ' 1 ' 

where c = arctan(/i)/ln(&), which exhibits the limit cycle. 

In the continuum limit of b — *■ 1, the sum over energy scales, X)i — * I dA/(Aln6). 
When the ratio r c = A/A of two cutoffs separated by one complete cycle is kept constant 
and b tends to 1, the difference p between powers and N must grow to infinity, as 
p = lnr c / In b. The corresponding coupling constant h = tan(7r/p) becomes then equal to 
7r/p = (-7T In b) I In r c and the constant c in Eq. (0) becomes equal to 7r/lnr c . This result 
shows that the limit b — > 1 in the discrete model reproduces the continuum case from 
Ref. @, and the period ratio is r c = exp(7r/c) for all values of the imaginary part c of 
the interaction term. The fixed-point behavior with asymptotic freedom is obtained only 
in the limit c — > 0, as the RG behavior in a one long cycle. 



3 Numerical example 



For illustration of the limit cycle from Fig. 1 and to show the existence of bound state 
solutions associated with it, Table 1 provides results of diagonalization of two Hamiltonian 
matrices: H\ with m and n ranging from M = —25 to iV = 16, and H 2 with N — 11 and 
the same M. Both Hamiltonians have the same values of b = 2, = hn = tan(7r/5) ~ 
0.726543, and giQ = gn = 0.060606 ||. The Hamiltonian H 2 is considered a result of five 
discrete RG steps applied to the Hamiltonian Hi. The infrared cutoff M was introduced 
in order to use a computer for calculations. The lower is M the better is the agreement 
between the example and the cycle theory applicable for M — > — 00. The eigenvalues 
of small modulus (not too small, however, to stay away from the lower bound b M ) are 
identical in both cases within the displayed numerical accuracy. For these eigenvalues, 
the RG analysis of the previous section is fully confirmed. For eigenvalues close to the 
artificial lower bound b M , the finite size of the matrix changes the pattern in a certain 
way, which is not relevant here. 

The scaling phenomenon of interest is described by showing the ratios of ratios of 
successive eigenvalues for one matrix to two characteristic numbers. The negative eigen- 
values, that correspond to bound states (E n = b n are positive for all values of n), appear 
in a geometric sequence with quotient 32, with accuracy approaching several ppm in the 
best cases. Thus, there is one bound state per cycle among the solutions. The positive 
eigenvalues also appear in a geometric sequence, but with a quotient almost equal 2 5 / 4 , 
which can be understood on the basis of dimensional analysis without inspecting details 
of the solution. Namely, after p reduction steps one eliminates p states from the dynamics 
and obtains the Hamiltonian H2 that is equivalent to b~ p Hi = Hi/32. Since in every 
cycle among the p states only p — 1 have positive eigenvalues and the cycle repeats itself 
indefinitely, the ratio of the successive p— 1 eigenvalues (if such a common quotient exists) 
must satisfy the condition r p ~ l = b p in order to scale properly from one cycle to another. 
Therefore, r = b p ^ p ~ 1 \ and in the example one obtains the factor 32 1//4 . Deviations from 
that rule within a single cycle cannot be understood so simply but the calculation shows 
a repeating sequence within a cycle with all elements very close to r. 

However, the most prominent feature of the example is the appearance of one negative 
eigenvalue in every cycle. It requires further studies to be fully understood. Here it is 
noted that when the period p tends to infinity with Hn — > 0, one recovers the case 
from Ref. ||, with asymptotic freedom and a bound state, but in infinitely many copies 
connected together in a sequence. In every single long-cycle link in that chain, looked at 
in the direction of growing A, the asymptotically free coupling decreases logarithmically 
towards large energies until it eventually crosses zero, grows in size and turns up rapidly. 
A new generation bound state is uncovered at that scale, and a new cycle starts from a 
large coupling constant that decreases slowly again. 
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n 


A = 2 16 


r 


A = 2 n 


r 


13 


0.954734+05 


- 


- 


- 


12 


0.328953+05 


0.819481 


- 


- 


11 


0.131845+05 


0.953270 


- 


- 


10 


0.545303+04 


0.983701 


- 


- 


9 


0.228087+04 


0.994830 


0.298354+04 


- 


8 


0.956198+03 


0.997093 


0.102798+04 


0.819481 


7 


0.401063+03 


0.997590 


0.412015+03 


0.953270 


6 


0.168593+03 


0.999803 


0.170407+03 


0.983701 


5 


0.709615+02 


1.001090 


0.712770+02 


0.994830 


4 


0.298253+02 


0.999652 


0.298812+02 


0.997093 


3 


0.125234+02 


0.998679 


0.125332+02 


0.997590 


2 


0.526682+01 


1.000260 


0.526853+01 


0.999803 


1 


0.221724+01 


1.001270 


0.221755+01 


1.001090 





0.931986+00 


0.999731 


0.932040+00 


0.999652 


-1 


0.391347+00 


0.998713 


0.391357+00 


0.998679 


-2 


0.164586+00 


1.000270 


0.164588+00 


1.000260 


-3 


0.692886-01 


1.001280 


0.692889-01 


1.001270 


-4 


0.291245-01 


0.999733 


0.291245-01 


0.999731 


-5 


0.122296-01 


0.998714 


0.122296-01 


0.998713 


-6 


0.514332-02 


1.000270 


0.514332-02 


1.000270 


-7 


0.216526-02 


1.001280 


0.216526-02 


1.001280 


-8 


0.910135-03 


0.999730 


0.910135-03 


0.999730 


-9 


0.382169-03 


0.998705 


0.382169-03 


0.998705 


-10 


0.160723-03 


1.000250 


0.160723-03 


1.000250 


-11 


0.676587-04 


1.001230 


0.676587-04 


1.001230 


-12 


0.284359-04 


0.999610 


0.284359-04 


0.999610 


-13 


0.119370-04 


0.998423 


0.119370-04 


0.998423 


-14 


0.501683-05 


0.999594 


0.501683-05 


0.999594 


-15 


0.210854-05 


0.999630 


0.210854-05 


0.999630 


-16 


0.882705-06 


0.995684 


0.882705-06 


0.995684 


-17 


0.367003-06 


0.988875 


0.367003-06 


0.988875 


-18 


0.150481-06 


0.975210 


0.150481-06 


0.975210 


-19 


0.584264-07 


0.923457 


0.584264-07 


0.923457 


-22 


0.654169-08 


0.266298 


0.654169-08 


0.266298 


-4 


-.622117-06 


0.964026 


-.622117-06 


0.964026 


-3 


-.206506-04 


0.998881 


-.206506-04 


0.998880 


-2 


-.661561-03 


0.999965 


-.661561-03 


0.999960 


-1 


-.211707-01 


0.999994 


-.211708-01 


0.999831 





-.677466-00 


0.999832 


-.677580+00 


0.994617 


1 


-.216826+02 


0.994617 


-.217999+02 


0.824124 


2 


-.697598+03 


0.824124 


-.846472+03 




3 


-.270871+05 









Table 1. The columns contain: n - the approximate integer powers of 2 5 / 4 (2 5 ) that give positive 
(negative) eigenvalues; A = 2 16 - the eigenvalues of the Hamiltonian with N — 16 and M = —25 (see 
the text for details); r - the ratio of the ratio of two successive eigenvalues to 2 5 / 4 (2 5 ) for the positive 
(negative) eigenvalues; A = 2 11 - the eigenvalues of H after 5 discrete RG steps (the same M is used). 
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4 Conclusion 



The possibility of a long asymptotically free cycle due to small imaginary couplings, with 
an abrupt end and a bound state formation where a new cycle starts, is an interesting 
feature of our specific model. But the central result of our paper is that such a simple 
Hamiltonian exhibits limit cycle behavior for non-zero h instead of a fixed point. Together 
with the limit cycle found by Bedaque et al, and implicit in earlier work of Efimov and 
Thomas, this raises the question of whether limit cycle behavior is a common feature 
of renormalized quantum theories, outside of those that have already been solved and 
are known to exhibit fixed points. The model analysis carried out here certainly needs 
improvements, especially concerning details of the bound state formation and the sharp 
turn of the coupling where it approaches 1 and Eq. (|) produces a jump. 

Separately from the existence of the limit cycle, the model is also valuable because it 
can be studied using the similarity RG technique (see [§]), and the latter can be applied to 
Hamiltonians in quantum field theories of major interest, such as QCD (see e.g. |], |IO|j). 
(We do not expect QCD in particular to show limit cycle behavior: our comments about 
limit cycles apply more to broad classes of Hamiltonians, most not yet examined, rather 
than specific ones already examined as extensively as QCD has been.) The authors have 
generated numerical solutions to Wegner's differential equation in the model and found 
continuous similarity RG evolution of effective couplings with all key features described 
here. Thus, when one uses the similarity scheme in studies of effective Hamiltonians in 
more advanced models, one may also encounter complex RG scenarios. It remains to be 
seen if the model structure reappears in realistic quantum theories. 
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